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Abstract 

At the NASA Lewis Research Center's Icing 
Research Tunnel ice shapes, similar to those 
which develop in-flight icing conditions, were 
formed on an airfoil. Under cold room condi- 
tions these experimental samples were care- 
fully removed from the airfoil, sliced into thin 
sections, and x-rayed. The resulting microra- 
diographs were developed and the film digit- 
ized using a high resolution scanner to extract 
fine detail in the radiographs. A procedure 
was devised to calibrate the scanner and to 
maintain repeatability during the experiment. 

The techniques of image acquisition and 
analysis provide accurate local density meas- 
urements and reveal the internal characteristics 
of the accreted ice with greater detail. This 
paper will discuss the methodology by which 
these samples were prepared with emphasis on 
the digital imaging techniques. 

Introduction 

Part of the Icing Research Program at NASA 
Lewis Research Center involves the study and 
measurement of local and total densities of ice 
shapes accreted on airfoils. Density meas- 
urements are needed to study the structure and 
formation of ice and to evaluate density corre- 
lations employed in computer applications that 
model ice accretion. 

This report presents the imaging methodology 
employed to obtain local densities in ice 
shapes accreted on airfoils at the NASA Lewis 
Icing Research Tunnel (IRT). The methodol- 
ogy takes advantage of modem imaging tech- 
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niques to improve and extend F. Prodi's 1,2 x- 
ray contact microradiography method of 
measuring ice density. In the past, measure- 
ments of accreted ice density have mainly been 
performed using the rotating multi-cylinder 
technique, 3,4 but there are difficulties in ex- 
tending the results to non-rotating bodies such 
as airfoils 5 . A methodology is needed to 
measure local and total densities of accreted 
ice shapes on airfoils. Prodi's technique, al- 
though originally applied to hailstones, can be 
improved and extended to measure local and 
total densities in ice shapes accreted on airfoils 
by using recent technical advances in imaging. 

The methodology described in this paper was 
applied to ice shapes that were accreted on a 
NACA 0012 airfoil in the IRT for a range of 
icing conditions. The shapes were then re- 
moved from the airfoil and cut into thin slices. 
X-rays were then used to obtain a contact mi- 
croradiograph of each slice. The developed 
film was digitized and image processing tech- 
niques were used to measure the density of the 
film which was correlated to the physical 
density of the ice (see Figure 1). 

Applications of the methodology are presented 
for a rime ice, glaze ice, and a mixed ice case. 
For the mixed and rime ice cases local densi- 
ties are measured across a region of the ice 
shape and image processing techniques are 
used to obtain information on the structure of 
the ice. 


Experimental 

Ice shapes were accreted on an NACA 0012 
airfoil in the Icing Research Tunnel 6 (see 
Figure 2) at the NASA Lewis Research Center 
for a range of icing conditions. The ice shapes 
were removed from the airfoil, wrapped in 
plastic and stored in a freezer. When enough 
ice samples were collected, the temperature in 
the tunnel was lowered to 0° F, the air flow 
stopped and the tunnel used as a cold room. 
The experimental hardware was moved inside 
the tunnel and allowed to chill until it reached 
tunnel ambient temperature. The ice samples 
were then removed from the freezer and cut to 
approximately 10 mm thickness with a band 
saw. 



Figure 1: Experimental schematic 


parallel. As the sample neared the 3 mm depth 
it was measured at eight locations around the 
perimeter with a depth gage until the desired 
thickness was obtained. 

Once the sample was determined to be 3 mm it 
was removed from the mounting block. The 
instant the sample became free from the block 
it was carefully placed on a sheet of mylar 
used to transfer it to the x-ray source. The x- 
ray device was fabricated from a photographic 
enlarger with an attached easel base (see 
Figure 3). A Watkin-Johnson x-ray source 
was used with a maximum output power of 50 
kV and a maximum output current of 1 mA 
was used. For safety reasons the x-ray device 
was shielded on three sides with 1/4 inch 
aluminum plates to block scattered rays. 





To mount the thick slices onto the microtome 
a method was developed in which an alumi- 
num plate was placed on a heating pad and an 
icy fingerprint was imprinted on it. When the 
plate reached the temperature where the fin- 
gerprint was beginning to melt the plate was 
removed and the sample placed on it. The 
plate was immediately placed on a large steel 
block which acted as a cold sink. The thin 
liquid film froze immediately creating a very 
strong bond. Extreme care was exercised to 
keep the thin film to a minimum to avoid con- 
tamination of the sample. 

One side of the sample was microtomed creat- 
ing a flat surface, then the sample was flipped 
over and the reverse side microtomed. Since 
the mounting plate was parallel to the cutting 
blade both surfaces of the ice sample were also 


Figure 2: Schematic of Icing Research Tunnel 

Kodak Industrex ready pack II type M x-ray 
film was used because of its high resolution 
and sealed packaging. The film was held in 
place by the easel and the sample placed in 
contact with the film packet. Prior to expo- 
sure the x-ray device was run for several min- 
utes to warm up and reach a stable operational 
mode. Exposure was determined in the pre- 
test phase of the experiment ensuring the 
density of the film was within the dynamic 
range of the film scanner. An input current of 
1 mA and 25 kV exposing the film for 30 sec- 
onds produced sufficient density. 

A reference sample of clear ice made from 
distilled, de-ionized water was placed in the 


comer of each radiograph exposed. This ref- 
erence sample was also microtomed to 3 mm 
with co-planar sides and served to calibrate 
the digital images. 



Figure 3: X-ray device 


Non-uniformities in the x-ray source caused 
density variations in the film which would af- 
fect the localized density measurements. To 
compensate a series of flat field images (see 
Figure 4) were exposed each day of the ex- 
periment. The flat field image is exposed under 
the x-ray source without a sample and digitally 
subtracted from the sample image (see Figure 
5). During the experiment, flat fields were ex- 
posed in between each sample image to 
monitor shifts in the non-uniformity and to 
provide an average flat field over time. Proc- 
essing of the film was performed with a Kodak 
X-Omat radiographic film processor at the 
recommended time and temperature. 



Figure 4: Flat field image 


The radiographs were scanned with a 
Leafscan-45 film scanner at a resolution of 
1 100 pixels per inch (ppi) creating a 19 MB 
file. Four v-shaped registration marks were 
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exposed on each sheet of film facilitating reg- 
istration of each image. A registration mark 
on each side of the film served as a backup if 
one shifted. By having several registration 
marks any two or all could be used to insure 
registration accuracy. 

Prior to the experiment the film scanner was 
tested to evaluate its consistency from scan to 
scan. It was found that the scanner auto- 
ranged the output gray values based on the 
radiograph density range. The auto-ranging 
changed the internal lookup table based on the 
minimum and maximum intensity value of the 
pre-scanned image, thereby expanding the dy- 
namic range of the output image. To circum- 
vent the auto-ranging feature a balanced target 
was pre-scanned, in this case a Kodak gray 
scale, which provided the same minimum and 
maximum intensity value for each subsequent 
scan. 



Figure 5: Sample Image with subtracted flat field 


Registration was performed using two of the 
four registration marks (see Figure 6) and 
checked against the other two. Each registra- 
tion mark was enlarged sufficiently to locate 
the end point. This was determined accurately 
by counting the cross-sectional profile of pix- 
els across the registration point and selecting 
the midpoint pixel. Then counting the pixels 
from the edge of the film to the same midpoint 
pixel and using these values as a standard 
measurement for the entire experiment. 
Shifting the images so each midpoint pixel re- 
sided at the same spatial location effectively 
registered each image to another. The dis- 
placement of opposite registration marks were 
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checked to verify proper alignment when 
scanned. If this varied by more than five pix- 
els the radiograph was realigned in the film 
holder and rescanned. Other point to point 
calculations were used to monitor shrinkage 
and swelling of the film which was found to be 
negligible. 



Figure 6: Magnified registration mark 


ice sample. Several hundred pixels from the 
reference sample region were averaged result- 
ing in a single intensity value. The conversion 
process was accomplished by making a global 
adjustment to the digital image, setting the 
average intensity value b of the reference 
sample equal to 0.917 gm cm 3 resulting in the 
corrected density image p(x,y). 


p(x,y) = 0.917 


(*.y) 

b 


(3) 


where b is the average intensity value of the 
reference sample represented by 




(4) 


Once all of the radiographs were scanned and 
registered, the flat field images for that par- 
ticular day were averaged together into a sin- 
gle flat field image g(x,y). Where fd t (x,y ) 

is each flat field image. 


This procedure converted all of the raw in- 
tensity values of the digital image into values 
that represent the actual ice sample densities. 

Results 


g(x,y) = -^^fd i (x,y) (1) 

The averaged flat field image was then sub- 
tracted from all of the sample images 5(x,y) 
which resulted in the sample image s' (x, y) 
with corrected intensity values. 

s’(x,y) = s(x,y)-g(x,y) (2) 

To correlate the system from physical ice 
density to digital intensity values there must 
exist a known density region in the x-rayed 
image. Each sample was exposed with the 
reference ice sample for calibration. Since the 
reference sample was produced with extremely 
clean water it was assigned the value of 0.917 
gm cm' 3 . This assumption served as the corre- 
lation point for all of the digital images in the 
experiment. 


The microradiographs reveal density variations 
within the sample that would have normally 
been concealed from the naked eye. Digitizing 
the samples allowed image processing tech- 
niques to evaluate ice formation and morpho- 
logical anomalies. Since the density of each 
pixel, as it relates to the reference sample, was 
known the image could be thresholded at 
specific densities revealing unobserved charac- 
teristics. Figure 7 shows a sample of glaze ice 
containing air bubbles. 



Figure 7: Image of glaze ice sample 


The density value assigned to the reference 
sample made it possible to associate the digital 
intensity values to the physical density of the 
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Figure 8: Thresholded sample to show air pockets 
( glaze ice) 



Figure 9: Glaze ice image with edge detection 

Figure 8 shows the same image thresholded to 
show the morphology of the sample. The ad- 
vantage of this technique allows one to ana- 
lyze the images with a resolution up to 256 
distinct densities and a spatial resolution of 
1 1 00 density measurements per inch. Figure 9 
shows how edge detection can be utilized to 
evaluate air pocket characteristics. 



Figure 10: Image of mixed ice sample 

Figure 10 is a typical example of a mixed ice 
sample. The thresholded mixed ice sample in 
Figure 1 1 shows air pockets along the periph- 
ery. 



Figure 11: Thresholded sample to show air pockets 
(mixed ice) 



Figure 12: Image of rime ice sample with profile line 

The rime ice sample in Figure 12 shows dis- 
tinct morphological characteristics. The in- 
tensity profile in Figure 13 shows that these 
are portions with lower density values and in- 
dicate air bubbles in the ice structure. Figure 
14 shows the thresholded rime ice sample. 



Figure 13: Profile of density values for rime ice 
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DISCUSSION 



Figure 14: Thresholded rime ice sample 


This experiment demonstrates that x-rays and 
digital image processing can improve the 
analysis of density variations in accreted ice 
samples. Using the method described it was 
possible to analyze the density variations 
throughout the sample and identify character- 
istics imperceptible prior to this experiment. 
The improvement in resolution over Prodi's 
method significantly enhances the data gather- 
ing capabilities. Prodi used a microdensitome- 
ter with an aperture size of 0.4 mm and a 
sampling interval of 0.5 mm. In this experi- 
ment each pixel represented a 23 micron aper- 
ture size with an evenly spaced 23 micron 
sampling interval. Our spatial resolution was 
17 times finer than the microdensitometer used 
by Prodi. 

In addition to improving resolution, digital 
imaging provides repeatability in measure- 
ments and absolute positional accuracy. With 
this ability verifying a measurement or analysis 
is simple and effective. Another significant 
aspect of digital imaging is eliminating the ne- 
cessity of handling the original radiographs for 
each measurement. The radiograph is scanned 
once and stored, leaving the digital image for 
analysis. All of the digital images were re- 
corded onto CD-ROM. This has the advan- 
tage of providing a convenient means to re- 
trieve, analyze, and archive the data. CD- 
ROM is advantageous for archiving because it 
bypasses the impermanence issues associated 
with silver halide systems. 


A digital method was applied to evaluate ac- 
creted ice densities by making local and aver- 
age densities measurements. This technique 
shows that the resolution of the imaging sys- 
tem is adequate to evaluate the structure and 
density characteristics of icing growth. It is 
felt that the use of digital imaging provides the 
ability to make repeatable and precise meas- 
urements. 

Variability can be attributed to basically two 
areas; the variation in thickness of the samples 
and the quantization of density values. The 
thickness of the sample was measured in eight 
locations around the perimeter during the mi- 
crotoming process. However, it was noticed 
that blade of the microtome flexed a minuscule 
amount from the mounted portion to the end 
of the blade when attempting to shave more 
than 1 5 microns in a single pass. This resulted 
in the sample having a sloped surface. By de- 
creasing the amount of distance the blade 
dropped, thus shaving less of the sample on 
each pass, the blade flex was reduced almost 
to the point of elimination providing parallel 
surfaces. 

Quantization error could effectively be re- 
duced by utilizing more than 8 bits of infor- 
mation, thus providing a wider dynamic range. 
The Leafscan-45 film scanner outputs 8 and 
16 bits of information which results in a dy- 
namic range of 256, and 65,536 density values 
respectively. 

Possible improvements to the experiment may 
be realized by utilizing an additional reference 
sample impervious to variability. This refer- 
ence sample, made from thin metal or plastic, 
would be used to monitor potential thickness 
variation in the ice reference sample. 
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APPENDIX 


Correlation of the density (blackness) of the film with the actual physical density of the ice sample is 
done following Prodi's methodology 2 . 

If we have a parallel monochromatic X-ray beam of intensity I o incident on the ice sample and on the 
film plate (Figure Al), we call X A the thickness of a region of the ice free of voids (air bubbles), p, the 
density of the ice in this region, I A the transmitted intensity, and D A the density of the film due to the 
exposure to the transmitted intensity. 

In a region of the ice containing air bubbles, X A is still the thickness of the ice sample (sample 
assumed to have parallel faces), X B is the true thickness of the material (thickness of the material less 
voids), p H is the density of the ice in this region, I B is the transmitted intensity, and D B is the density 
of the film due to the exposure to the transmitted intensity. 
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The ratio of the density can be related to the ratio of the thickness of the two regions: 

P„=ni B /X A (Al) 

where nig is the mass contained in volume X A . 1 . 1 in the region with air bubbles. But ir^ can be 
expressed as: 

®B = P. X B 

then p H = p, X B /X A 

and p H /p, = X B /X A (A2) 

The ratio of the thickness X A / X B can be related to the beam intensity using Lambert's law: 

I = I 0 exp ( - G X ) (A3) 
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where I o and I are the incident and transmitted radiation respectively, a is the absorption coefficient of 
the material, and X is the thickness. 

From (A3) X = -l/oln(I/I 0 ) (A4) 

and X B /X A = ln(I o /I B )/ln(I o /I A ) (A5) 

where the coefficient of absorption has been canceled out, assuming the ice is homogeneous wherever 
there are no bubbles. 

Substitution of (A5) into (A2) gives: 

p H / p, = ln(I 0 /I B )/ln(I 0 /I A ) (A6) 

Photographic exposure E can be expressed as: 


E = It (A7) 

where t is the exposure time. The response characteristic of the film for a certain photographic 
emulsion is defined as: 

y = D / log E = D / log (I t) = D / ( log I + log t ) (A8) 

where log I = ( D - y log t ) / y (A9) 

If we operate in the linear region of the film where y is constant, then 

p„ / p, = [((D 0 - y log t) / y ) - «D b - y log t) / y )] / [((D 0 - y log t) / y ) - ((D A - log t) / y )] 

p H / Pl = (D-D b )/(D 0 -D a ) 

where D o ,D A , D B , are respectively the measured photographic density in the region clear of the 
sample, in the region where the ice is clear of bubbles, and in the region where the ice contains 
bubbles. 

The local density of the ice sample in a given location can be obtained if we know the density at a 
location within the sample where the ice is clear of bubbles, and if we measure the film density at 
three locations: in a region where the ice is clear of bubbles, in a region of the film away from the 
sample, and in the region where we want to measure the physical density of the ice: 

Ph = P, (D - D b ) / (D o - DJ 
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